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Abstract 

We propose to introduce a new stochastic process in molecular dynamics in 

order to improve the description of the nucleon emission process from a hot 

nucleus. We give momentum fluctuations originating from the momentum 

width of the nucleon wave packet to the nucleon stochastically when it is 

being emitted from the nucleus. We show by calculating the liquid gas phase 

equilibrium in the case of antisymmetrized molecular dynamics, that with this 

improvement, we can recover the quantum mechanical statistical property of 

the nucleus for the particle emission process. 
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In multifragmentation reactions in intermediate energy heavy ion collisions, many inter- 
mediate mass fragments are produced. Data of the mass distribution and other exclusive 
quantities have been reproduced rather well by the statistical models which assume 
the decay of the thermal source into fragments according to the statistical weights of the 
possible final states. It is therefore reasonably considered that the multifragmentation pro- 
cess is largely governed by the statistical property. On the other hand, molecular dynamics 
simulations try to describe the fragmentation based on the microscopic dynamical evolution 
of the system |3|-f|. With molecular dynamics simulations we can study the dynamical and 
statistical aspects in a unified framework without assuming any thermal equilibrium. How- 
ever, for the reliable description of the fragmentation, we should check the validity of the 
statistical aspect of molecular dynamics and, if necessary, modify the dynamics of molecular 
dynamics so as to get the appropriate statistical aspect. 

In our previous studies |J|7]] of heavy ion collisions with antisymmetrized molecu- 
lar dynamics (AMD) in the incident energy region between 30 MeV/nucleon and 100 
MeV/nucleon, we observed that light nuclei such as 12 C break up into fragments easily 
in the dynamical stage of the reaction but intermediate mass fragments are seldom pro- 
duced in the collision of heavier nuclei such as 40 Ar. We are also aware that the decay of 
excited nuclei produced in the reaction is too slow when the AMD calculation is continued 
for a long time, compared to the prediction by the statistical decay model. These points 
strongly suggested that we should improve particle emission processes in AMD. In this let- 
ter, we propose a method how to improve AMD in the nucleon emission process from the 
excited nucleus and show by calculating the liquid gas phase equilibrium that the improved 
AMD describes the quantum mechanical statistical property of hot nuclei for the particle 
emission process while the original AMD describes only the classical statistical property for 
that process. Our method of improvement is of general nature and can be applied to other 
molecular dynamics simulations. 

Before the discussion of the statistical aspect and the improvement of AMD, we will 
explain the usual AMD || very briefly for the convenience of the readers. AMD describes 
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the nuclear many body system by a Slater determinant of Gaussian wave packets as 

$(Z) = det[exp{-v( rj - Zi/y/U)^^)], (1) 

where the complex variables Z = {Zj} are the centers of the wave packets. We took the 
width parameter v = 0.16 fm~ 2 and the spin isospin states Xa t — P t> P I, n t> or n I- The 
time evolution of Z is determined by the time-dependent variational principle and the two- 
nucleon collision process. The equation of motion for Z derived from the time- dependent 
variational principle is 

Ci a j T with a, t = x, y, z is a hermitian matrix, and 7i is the expectation value of the 
Hamiltonian after the subtraction of the spurious kinetic energy of the zero-point oscillation 
of the center-of-masses of fragments, 

n = (H)-^A + T (A-N F ), (3) 

where Np is the fragment number and T is 3hu/2M in principle, but treated as a free 
parameter for the adjustment of the binding energies. In two-nucleon collisions, physical 
coordinates Wj are introduced as 

Wi-^Q)^ Qij - dm . z . } — • (4) 

Except for the two-nucleon collision process, the AMD equation of motion (0) constitutes 
a kind of Hamilton system, and therefore one may think that the statistical property of 
AMD is not quantum mechanical ||. It is not true because AMD can give the exact time 
evolution of the wave function for the harmonic oscillator mean field. Namely although the 
motion of the wave packet centers is classical, the AMD wave function containing the spatial 
and momentum spread of wave packets can have quantum mechanical statistical property. 
However, for the process of particle emission, the statistical property of AMD may be of 
classical nature, because the AMD description of the particle emission does not duely take 



into account the momentum spread of wave packets. In order to get physical picture and to 
enable further development, we calculate the statistical property of the excited nucleus for 
the particle emission process in a new method explained in the following. We added to the 
Hamiltonian (|3]) the potential wall with a large radius of the form 

~£/(|D;-D C m|) (5) 

z i 

with 

f(x) — (x — a) 2 9(x — a), 

y/i>T>i = ReZi, D C m = — EjDj, 

•Atot 

a = 12fm, fc = 5MeV/fm 2 , 

and put A tot nucleons into the potential wall and give the total energy E tot . The time 
evolution is calculated for a long time (~ 10 5 fm/c). As shown in Fig. [l| there are usually a 
nucleus and nucleons in the potential wall. Nucleons are sometimes emitted and absorbed 
by the nucleus. This situation can be interpreted as the phase equilibrium of liquid (nucleus) 
and gas (nucleons). Typical values of the mass number of the largest fragment An q and its 
internal energy E\ iq depend on the choice of A tot and -Etot- We select the moments at which 
An q takes a given value in order to get the statistical property of the nucleus of the given 
mass number, and the long time average value of Eu q and the temperature are calculated. 
The temperature T, which should be common to both phases, is calculated as the long time 
average value of r, where (3/2)r is the kinetic energy (plus the potential energy from the 
wall) per nucleon in gas phase. Since the gas phase is dilute, the effect of the Pauli principle 
is neglected. When there are small fragments in the gas phase, they are excluded in the 
calculation of the temperature. The calculated relation between E\ iq and T is shown in the 
left part of Fig. |2] by squares, which lie just on the line {E\ iq )/A\ iq = 3T — B.E. of the classical 
statistics. Each calculated point corresponds to a specific choice of Atot and -Etot- In the 
calculation of the temperature T, i.e., the kinetic energy of gas nucleons, we only used the 
central value of the wave packet Im Zj in momentum space. The neglect of the momentum 

4 



00 

FIG. 1. The density snapshots of the many-nucleon system contained in a potential wall 
(shown by circles) calculated with the usual AMD. These figures are for the total number of the 
nucleons Atot = 15 and the total energy Etot = — 50 MeV, and the calculated temperature and the 
energy of the nucleus are T = 1.48 MeV and {E^ q ) /A\i„ = —3.65 MeV for the mass number of the 
nucleus A\i q = 14. The Volkov force is used for this calculation. 

width is consistent with the subtraction of the spurious zero-point kinetic energy from the 
Hamiltonian (Q). If we consider the momentum width, the temperature will shift by (2/3)Xo 
(~ 6 MeV) and there will not be any states with temperature less than (2/3)T . 

Microscopic processes which are important for the liquid gas phase equilibrium is the 
emission and the absorption of nucleons by the nucleus. The above result indicates that 
AMD has a problem in the description of the nucleon emission from the excited nucleus, 
since nucleon absorption is naturally described by AMD. Each nucleon in the nucleus has 
a momentum width of the wave packet which is an important part of its Fermi motion in 
the nucleus. When such a nucleon is emitted from the nucleus in AMD, it will have the 
momentum corresponding to the central value of its momentum distribution in the nucleus. 
Especially if the central value is not sufficiently large, the nucleon will never be able to go 
out, even if there should be some probability of going out due to the tail of the momentum 
distribution. In order to cure this problem, we regard the emitted nucleon as a classical 
particle without momentum distribution or a plane wave, though it is a wave packet in a 
nucleus. And we add a new stochastic process of giving the momentum fluctuation of the 
wave packet to a nucleon when it is being emitted from a nucleus. The total momentum and 
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FIG. 2. The statistical property of the excited nucleus calculated with the usual AMD and the 

AMD-MF. The adopted effective interaction is the Volkov force in the left part, and the Gogny force 

in the right part. Lines of (E/A) hq = 3T + const., T 2 /(12 MeV) + const., and T 2 /(8 MeV) + const. 

are drawn for the comparison. 




FIG. 3. Similar to Fig. but calculated with the AMD-MF. A tot = 28 and E tot = 60 MeV 
for these figures, and calculated temperature and the energy of the nucleus are T = 6.12 MeV and 
(E\i q ) /A\i q = —5.52 MeV for the mass number of the nucleus ^4n q = 14. 
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energy is conserved by adjusting the state of the nucleus. This new version of AMD is called 
AMD-MF and the details will be described in following paragraphs. By this incorporation 
of the momentum fluctuation into AMD, the situation of the liquid-gas phase equilibrium 
changes as shown in Fig. |3|. Even when the excitation energy of the nucleus is smaller than 
in the case of usual AMD (Fig. [TJ) , the number and the averaged energy of the gas nucleons 
are much larger. The AMD-MF is not a Hamilton system any more, and the statistical 
property of the excited nucleus for the particle emission process has changed drastically into 
quantum mechanical one, as shown in Fig. § by diamonds, which is similar to the empirical 
Fermi gas formula (E hq )/A hq = T 2 /(12MeV) - B.E.. 

It should be noted that our statistical ensemble for an excited nucleus is between the 
canonical ensemble and the microcanonical ensemble. Although both of E\ iq and r are 
fluctuating in time, they are correlated due to the energy conservation relation 

AE hq + ^A gas Ar = (6) 

with A gas = A tot — A\i q , as long as one can neglect the probability of the existence of 
light fragments in gas phases and the small interaction energy between both phases and 
among the gas nucleons. By choosing a large potential wall so that A gas is very large, it 
is possible to get the canonical ensemble with small At. On the other hand, if one take 
small volume for gas phase, microcanonical ensemble with small AE\ iq can be obtained. 
With our choice of parameters for the case of Fig. |3], for example, the energy fluctuation is 
found to be ^[AE^)/A\ iq = 1.2 MeV, which corresponds to the fluctuation of temperature 
^/(Ar 2 ) = 0.8 MeV. 

In AMD-MF, the momentum fluctuation is given to a nucleon when it is being emitted 
from a nucleus. At each time step, it is judged whether each nucleon is in a nucleus or it is 
isolated in free space. For each nucleon i, a representative point r-j is randomly generated 
around ReWi/y^ according to the Gaussian distribution, where W, is the physical coor- 
dinate (^). The density at without self-contribution is calculated in an approximate way 
as 
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/2lA 3 / 2 r n 

exp[-2( v ^r i - ReW/j . (7) 

If pi < O.lpo with p = 0.17 fm~ 3 , the nucleoli is considered to be isolated in free space; 
else it is in a nucleus. In the same way, it is checked whether at the previous time step the 
representative point — Re(Wj(t) — Wj(£ — At)) was in a nucleus or isolated. Momentum 
fluctuation is given to the nucleon i if it was in a nucleus at the previous time step and it is 
now isolated in free space. 

When it is decided that the nucleon i should be given the momentum fluctuation now, 
the following procedures are made. First of all, the nucleus from which the nucleon i is 
being emitted is defined as the cluster judged by the chain clustering method. The mass 
number of the nucleus is denoted by A nuc , which do not include the nucleon i. Then the 
relative momentum between the nucleon % and the nucleus is decided stochastically, by taking 
account of the fact that the momentum distribution of the nucleon and the distribution of the 
momentum per nucleon of the center-of-mass of the nucleus are approximately represented 
by Gaussian distributions 



exp 



respectively, where 



p-P i ) 2 /2h 2 u and exp -A nuc (p-P nuc ) 2 /2^ , (8) 



P, = 2/i^ImW i and P„ uc = E ImW r ( 9 ) 

nuc j'gnuc 

In order to keep consistency with the method of subtraction of spurious center-of-mass 
motion of fragments (Q), the energy To is subtracted from the relative kinetic energy by 
reducing the relative momentum. If the relative momentum is too small for this procedure, 
this momentum fluctuation is canceled. 

The first candidate of the final state of the momentum fluctuation process is determined 
by changing only the relative physical momentum between the nucleon % and the nucleus 
to the value decided in the previous paragraph. However, since the energy is not conserved 
with this state, the energy is adjusted to the initial energy with the least modification of the 
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state by applying the frictional cooling/heating method with A + ifj, — for the limited 
coordinates Zj with j in the nucleus. And after each step of the friction, the physical 
coordinate of the nucleon i and the physical center-of-mass and the total momentum of the 
nucleus are restored to the values of the first candidate. 

There are several possibilities for the above procedure to fail. The friction leads to a 
converged energy but still may not reach the initial energy. Or the transformation from 
W variables to Z variables may not exist due to the Pauli-blocking. In these cases, the 
momentum fluctuation is tried again by generating a different sample of relative momentum. 

In summary, a new stochastic process has been added to AMD in order to improve the 
description of the nucleon emission process from a hot nucleus. The momentum fluctuation 
originating from the momentum width of the nucleon wave packet in the nucleus is given to 
the nucleon stochastically when it is being emitted from the nucleus (AMD-MF). With this 
improvement, the quantum mechanical statistical property of the excited nucleus for the 
particle emission process has been obtained, which are calculated in the phase equilibrium 
of the liquid phase and the gas phase of nucleons. It should be noted that we need not 
totally modify the dynamics of the AMD equation of motion. As we mentioned before, we 
can believe that the AMD wave function gives fairly good time-evolution of the hot nucleus, 
because AMD can give the exact time evolution of the wave function for the harmonic 
oscillator mean field. What we have incorporated in this work is a minor branching process 
which is brought about by the tail of the nucleon wave packet. It may be necessary to 
generalize the modification so as to improve the cluster emission process in addition to the 
nucleon emission. For the application to heavy ion collisions, other processes caused by the 
wave packet tail should be considered such as the nucleon transfer between the projectile 
and the target. The wave packet width in spatial coordinate space may also play some role. 
Such extensions have been formulated and calculations of nuclear collisions have shown that 
fragmentation processes are largely influenced by the introduction of the new stochastic 
process described here and its extensions. We will discuss these in other papers. 
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